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Abstract 


We derive a systematic procedure to compute Green functions for the massive Schwinger 
model via a perturbation expansion in the fermion mass. The known exact solution of the 
massless Schwinger model is used as a starting point. 

We compute the corrections to the Schwinger mass up to second order. 
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1 Introduction 


The massless Schwinger model is wellknown to be equivalent to a theory of a free massive 
boson with Schwinger mass fi g = — ([|l[]). This massive field is formed via the chiral anomaly 
and may be interpreted as a fermion-antifermion bound state (f|, ||, fTo|, [Pj, |l9|). Besides, the 
massless Schwinger model shows other nontrivial features like fermion condensate, instantons 

(i i-n, 


In the massless 


and nontrivial vacuum structure (9 vacuum) (0-0, m-E, ta-yi, 

Schwinger model physical quantities do not depend on the vacuum angle 9. 

The massive Schwinger model is different in some repects. First, it is no longer exactly 
solvable (|l3|] |l(J, ]l(J, [f|). The physical state - the massive boson - is no longer free, and its 
mass acquires corrections due to the interaction. Instanton-like gauge fields and a nontrivial 
vacuum structure persist to be present, and, in addition, physical quantities now depend on 
the vacuum angle 9. The fermion condensate, too, acquires corrections due to the fermion mass 


(HI). 

Here we show how to compute Green functions for the massive Schwinger model within the 
Euclidean path integral formalism, using a perturbation expansion in the fermion mass. (The 
existence and finiteness of the mass perturbation theory for the massive Schwinger model was 
proven in jI0). For this purpose we use the known exact solution of the massless Schwinger 
model as a starting point. From the perturbation expansion we explicitly calculate corrections 
to the Schwinger mass up to second orcer in the fermion mass m, for a general vacuum angle 9. 


2 Exact solution of the massless case 

The vacuum functional (and Green functions) of the massive Schwinger model may be inferred 
from n-point functions of the massless Schwinger model via an expansion in the fermion mass. 
Indeed, we may write for the Euclidean vacuum functional (k... instanton number) 


Z(m,9)= J2 e ike Z k {m) 


( 1 ) 


k=— oo 


where 


Z k (m) = N 




oo n n 

m 


= N / D*D*D0 k £ —r IT / dx^ix^ixi 


„ n! . ,. 

n =0 i=l 


• exp 


dx 


+ -fd k a 2 p k } 


( 2 ) 


= e flu d u /3 corresponding to Lorentz gauge). Therefore the perturbative computation of 
Z(m,9) is traced back to the computation of scalar VEVs ( (1L S{ x i))o, S(x) = , F(x)T(x) ) 
for the massless Schwinger model and some space time integrations. It is useful to rewrite the 
scalar densities in terms of chiral ones, S(x) = S + (x) + S-(x), S± = because then only 

a definite instanton sector k = n + — contributes to the VEV (nT=i S+(. x i) U]ZiS-(xj)) 0 .A 
general VEV may be computed exactly (see e.g. §-@, @, @), 

{S Hl (x l )---S Hn (xn))o= (^) exp[^(-) CTi ^'47T D^Xi-Xj) 

i<j 


( 3 ) 
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where oi = ±1 for H[ = ±, D, J0 is the massive scalar propagator, 


1 


Dy. 0 {x) =-—K 0 (no\x\), D w ip) = 


-1 


P 2 + P o 


2 ’ 


(4) 


(K 0 ... McDonald function) and S is the fermion condensate of the massless Schwinger model, 


S = (**}o = 


(5) 


( 7 ... Euler constant). The index 0 for /j,q indicates that it is the order zero result, the index 
0 for the VEVs means that they are with respect to the massless Schwinger model. From this 


( 6 ) 


Z(m,9) may be computed (see [] 2 l| for details, |24| for its physical implications), 

Z{m,0) = e Va(m ’ e \ 

a(m, 9) = m —2 cos 9 + m 2 ^—^ (E cos 29 + F) + o(m 3 ) 

(V ... space time volume) where 

1 


E = J d 2 xE{x) = J d 2 x(e~ 2K °^ - 1 ) = -8.9139 ■ 

F = j d 2 xF{x) = J d 2 x(e 2Ko ^o\x\) - 1 ) = 9.7384 • \ 




(7) 


and for F a subtraction of a free field singularity is necessary. 

In order to compute VEVs for the massive Schwinger model one has to insert the corre¬ 
sponding operators into the path integral (1), (2) and divide by the vacuum functional Z(m, 9): 


(0) m = 


1 


(^EtII / dxi'&(xi)'&(x i )) 0 


Z{m,9) 


„ n\ , 1 

n =0 i=l 


( 8 ) 


Via the normalization all volume factors cancel completely, as it certainly has to be (we will 
explicitly see this in the computations). For the computation of VEVs of scalar (5) and chiral 
(S±) densities formula (3) suffices and could be used e.g. for the computation of the fermion 
condensate (\h'I') m (this however may be derived at once from ( 6 ), see [ 21 ]). 

For the computation of the Schwinger mass additional Green functions of the massless 
Schwinger model are needed. It is wellknown that the vector current correlator is the free 
massive propagator in the massless Schwinger model: 


1 


tw'dy F)y 0 {x y ) 

Therefore, for a perturbative calculation we need the n-point functions 

{My 2 ) J u(yi) n S Hi (xi))o = 


(9) 


i=l 


'E\«r 

-" 2 "/ 


Dy 0 (yi -y 2 ) + 4n(J2(-) (Ti D ll0 (x i — 2/1)) 

i=l j =1 

. e 4n I2k<i(-')' Tk ‘ TlD n( x k-xi) 


( 10 ) 


which may be computed from the Euclidean path integral or from bosonization (the latter 
method is easier). 
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3 Computation of the Schwinger mass 


For the perturbative computation of the massive propagator we simply have to insert successive 
orders of equ. (10) into ( 8 ). The factor (^e M is the same for all orders, therefore we 
will ignore it in the sequel. 

There are two terms from (10) to be inserted in first order, namely Sh,H = ±, and four 
terms, Sh 1 Sh 2 i = ±, in second order. For the moment we consider the case where all Hj = +; 
the other contributions can be inferred from this one by a rearrangement of some signs. Up to 
second order, we find ( {J^{y 2 )Ju{yi))m = ^ w '^'e w '^T( 2 /i, y 2 ) ) 


1 r X S /" 

T(yi,y 2 ) = Z ( m Q y [D^yi - y 2 ) + m-VD^iyx - y 2 ) + 4:TTm- J dxD^x - y 2 )D IXo (x - y x )+ 


9 

m 

~ 2 T 



J dx 1 dx 2 [D w (y 1 - y 2 ) + 4 tt(D w (x 1 - y 2 ) + D w (x 2 - y 2 ))- 
■ [D , 0 (xi - Vl ) + D w (x 2 - yi ))]e 4 ^o(-i- 2 )l 


( 11 ) 


Inserting equ. (6) for the vacuum functional Z(m,9 ) and using the perturbation formula of 
second order, 

ao + aqx + a 2 x 2 + o(x 3 ) 

1 + b\x + b 2 x 2 + o(x 3 ) 

ao + a\x + a 2 x 2 — a^bix — aob 2 x 2 — a\b\x 2 + aob 2 x 2 + o(x 3 ) (12) 


we arrive at 


T(yi,y 2 ) = D w (y 1 — y 2 ) + 47 rm— 


dxD fl0 (x - y^D^x - yi)+ 


TTTs^ 1 / ^ 1 \ 2 r 

4vr / ^idx 2 [2T> w (xi - y 2 )D w (x i - yi) + 2T> w (xi - yi)D^(x 2 - y 2 )\■ 
■{E{xi - x 2 ) + 1) - Anm 2 ^—) V / dxD fl0 (x - y 2 )D flQ (x - y x ) 

= D w (;yi -y 2 ) + 4irm— / dxD^ (x - y^D^ (x - yi)+ 


^ ^ 2 r 

m2 (^) 4:7tE J dxD w( x - V 2 )D w {x - yi)+ 

m 2 (—) 47 t / dxidx 2 D fJ , 0 (xi - yi)D^ 0 (x 2 - y 2 )(E{x 1 - x 2 ) + 1) (13) 

where E,E(x) are given in (7) and we used the x —> — x symmetry of all occurring functions. 
Observe that, as claimed, all volume factors V have dropped out. 

However, expression (13) is not yet the desired result for further computation, there is still 
one unwanted term. The last line of expression (13) consists of two terms, due to the factor 
(E(x i — x 2 ) + 1). The second one, proportional to 1, reads 


m 2 (^) 4vr J dxiD /X0 (x l -y l ) J dx 2 D fl0 (x 2 - y 2 ) 


(14) 
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and obviously does not contribute to the y\ — y 2 correlation. It is a disconnected part stemming 
from the VEVs 

(MVi) S +( x j))o ~ ^n'dy'iHuoixj - yi) (15) 

and must be subtracted from (13). 

Looking for the moment at the first order correction only, it is very easy to find the Schwinger 
mass correction: 

D yi tU) (Vi,V 2 ) = S(yi -y 2 ) + Mo^Vo (Vi ~ 2 / 2 )+ 

S f 

47 rm-(L> M0 (yi - 2 / 2 ) + Mo / dxD^x - y 2 )D^(x - yi)) 


where 


= <%i - 2 / 2 ) + Mi An (y\ - y 2 ) + o(m 2 ) 


p i = Mo + 47r?n( 


To obtain the second order result, we rewrite expression (13), without the disconnected part, 
in momentum space and substitute all functions by their Fourier transforms (thereby the con¬ 
volutions turn into products): 


— — 1 El 2 

T(p) = —z 9 + 47tto— —^ owT + 47rm 

P 2 + Mo 2 (p 2 + p 2 ) 2 


'E \2 1 


-i n , E 1 

- 9 (1 - 47 rm—- 9 

P 2 + Mo 2 p 2 + p 2 


TT—(E + E(p)) = 


2) {p 2 + pl) 




---=;-=---=;--- 1 - o(m 3 ) (18) 

p 2 + Mo + 47rm§ + 4vrm 2 (f ) 2 {E + E(p)) + (47rm| ) 2 pr^i 

Therefore, for finding the mass pole, p 2 has to obey the self consistency equation (after a 
rescaling p' 2 = E' = E{p q = 1) = p^E etc.) 


,2 , , m S 

p = —1 — 47T- 


m E , m 2 / E \2 47r , 

-47t —tt (- \E' + E'(p') + —5- 

/Jo 2/J.0 Mo ' 2/Jo 2 n' 4- 1 


p' +1 


The second order part (the term in square brackets) may be rewritten like 

[■■•] = | d 2 z[e- 2Ko(N) - 1 + e ip ' x (e- 2Ko{M) - 1 + 2 K 0 (\x\)] = 

n oo r2ir 

/ drr[2vr(e- 2Ko(r) - 1) + / kcos 0 ( e - 2 ^ o (r) - 1 + 2/i 0 (r))] = 

Vo Vo 

r OO 

2tt / drr[e" 2Ko(r) - 1 + J 0 (|p , |r)(e _21fo(r) - 1 + 2/i 0 (r))] (20) 

Jo 

where Jo is the Bessel function of the first kind. This expression behaves well around \p'\ = i 
and therefore we may set \p'\ = i because deviations from this value are of higher order in m. 
Using J 0 (r) = Jo(rr) we find 


,2 . m E 9 m” / i 

p' =-1-4 tt-87t 2 — 9 - (- • A 

Un 2Jin US V2iin ' 


m 2 / E \2 


Mo 2 Mo Mo V 2 Mo- 
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roc 

A := / drr[e~ 2K °^ — 1 + /o(r)(e~ 2A °^ — 1 + 2A"o(?’))] 
io 

A = -0.6599 (22) 

This result was computed for the positive chirality densities S+. For the pure negative chirality 
densities S- the result is completely identical. However in second order there are mixed products 
S + S_, S_S + , too. Both of them lead to the same result, where the above integral A is substituted 
by a similar expression B , 

r OO 

B:= drr[e +2K °^ — 1 + Io{r)(—e +2K °^ + 1 + 2it'o(?’))] 

Jo 


B = 1.7277 (23) 

In this expression the nice feature of cancellation of UV divergencies occurs. Indeed, both e 2K °^ 
and —Io(r)e 2K °( r ' 1 diverge like \ for small r (this divergency corresponds to the free fermion 
field divergency of the underlying theory), but obviously the divergencies cancel each other. 


In fact, this cancellation was already observed twenty years ago in [ 131 within a bosonization 
approach, and on general grounds it should continue to hold for higher orders. 

Collecting all results and multiplying each contribution with its corresponding 6 factor 


{S± 


0 ±id 


see (1-3)), we find for the Schwinger mass in second order 


— p 2 = rz A = 1 + 87t 

Mo 


— — cos 9 + 167r 2 ^ " {A cos 29 + B) 

Mo 2/xo Mo v 2 Mo ' 


(24) 


or, inserting all numbers (remember equ. (5)) 

m| = Mo(^ + 3.562! 


2 2 

771 TY1 777 

— cos e + 5.4807 ■ - 2.0933 • 

Mo Mo Mo 


cos 29) 


(25) 


which is our final result. 

For the special case 0 = 0 our result (24) precisely coincides with the result in [ 251, where 
the second order correction for 0 = 0 was computed within bosonization and using near light 
cone coordinates. In the same article this result was compared with a lattice calculation (|2G|]), 
and a good agreement is obtained within the range of the expansion parameter ^ where the 
lattice calculations were performed. 


4 Summary 

We have demonstrated a general method of computing n-point functions for the massive 
Schwinger model in mass perturbation theory within the Euclidean path integral formalism, 
using the known exact solution of the massless Schwinger model as a starting point. 

All n-point functions exist perturbatively and are finite in the infinite volume limit. 

Using this approach we were able to compute the mass perturbation corrections to the 
Schwinger mass up to second order, and almost all calculations could be done analytically. This 
feature does not persist to hold for higher orders. Already in third order an equation analogous 
to (19) will be a true self consistency equation that can be solved only numerically. 

For 0 = 0, the result (24) is well reproduced by lattice calculations ([f25|). 

The numerical calculations in this article were done with Mathematica 2.2. 
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